Time-dependent single electron tunneling through a shuttling nano-island 
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We offer a general approach to calculation of single-electron tunneling spectra and conductance of 
a shuttle oscillating between two half-metallic leads with fully spin polarized carriers. In this case the 
spin-flip processes are completely suppressed and the problem may be solved by means of canonical 
transformation, where the adiabatic component of the tunnel transparency is found exactly, whereas 
the non-adiabatic corrections can be taken into account perturbatively. Time-dependent corrections 
to the tunnel conductance of moving shuttle become noticeable at finite bias in the vicinity of the 
even/odd occupation boundary at the Coulomb diamond diagram. 
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I. INTRODUCTION 

Single electron tunneling (SET) is a salient feature of 
quantum transport in nanostructures. The SET phe- 
nomenon is observed in various systems, e.g. quan- 
tum dots in a tunnel contact with metallic electrodes?^ 
molecular bridges between the edges of broken metallic 



wires 
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atoms and molecules absorbed on metallic sur- 
faces in a contact with the tip of tunnel microscope^ 
etc. The study of electron tunneling through the nano- 
object with time-dependent characteristics is one of the 
most challenging problems in this field. 

There are several sources of time dependence, which 
may be realized in practical devices. The simplest one 
is the time-dependent gate voltage v g {t) applied to the 
dot. It is well known&2*i£ that this time dependence 
may be converted into the time dependence of tunnel ma- 
trix element. Another possibility is the nanoelectrome- 
chanical shuttling (NEMS)jii where the nano-size island 
suspended on a pivot— or attached to a string^ oscil- 
lates between the leads under the action of an electro- 
mechanical force. In case of molecular bridges, the vi- 
bration eigen modes may be the source of the periodical 
oscillations of tunneling parameters 

Usually the tunneling between metallic leads and such 
a nanoobject is accompanied by many-particle Kondo 
screening effec t 15 ! 16 resulting in specific type of zero-bias 
anomaly (ZBA) in tunnel conductance. Modification of 
Kondo regime because of periodically modulated in time 
tunneling rate due to the center-of-mass oscillations was 
studied recently in several papers. If the oscillations are 
the eigen modes of a nanoobject (molecule), then the 
Kondo-peak (zero bias anomaly in tunnel conductance) 
may transform into dip due to the destructive interfer- 
ence with vibrational mode— In case of adiabatic mo- 
tion induced by electromechanical forces (NEMS)ii the 
Kondo temperature follows the periodical evolution of 
the dot position and increases eventually due to effective 
reduction of the average distance between the dot and 
the leads,— which is determined by the mean square dis- 
placement of the dot (in analogy with Debye- Waller ef- 
fect in scattering intensity) . Non-adiabatic enhancement 



of Kondo tunneling through such moving nanoobject at 
finite source-drain bias has been studied recently^ 8 

In the present paper we consider adiabatic and non- 
adiabatic time-dependent effects in conventional cotun- 
neling between metallic leads due to periodic modulation 
of lead-dot tunneling rate. To suppress many-particle 
Kondo screening effects, one should consider leads with 
magnetically polarized electrons. The tunneling be- 
tween the ferromagnetically ordered leads was discussed 
recentlyi^^i^i^ 3 .^^, in the context of Kondo effect. 
We are interested in the situation, where the spin-flip co- 
tunneling is completely suppressed at small lead-dot bias 
and low temperatures. Such tunneling regime may be 
realized in half-metallic ferromagnets, where the Fermi 
surface is formed only by the majority spin electrons, 
whereas the spectrum of minority spin carries is gapped. 
The electronic and magnetic pro pert ies of such metallic 
compounds are reviewed in Ref. |2q . From the point of 
view of existing devices, where the leads are formed by 
two-dimensional electrons in degenerate semiconductors, 
the relevant material for our studies is dilute magnetic 
semiconductor (Ga,Mn)As— The indirect magnetic ex- 
change between Mn impurity ions is responsible for the 
long-range ferromagnetic order in this material. This ex- 
change is mediated by spin-polarized carriers near the 
top of the valence band. The tunnel current in the half- 
metallic regime arises due to the minority spin hole co- 
tunneling. 

We will show that in the absence of Kondo effect the 
problem of tunneling through moving nano-island (quan- 
tum dot) may be solved by means of time- dependent 
canonical transformation, which exactly takes into ac- 
count both adiabatic and non-adiabatic lead-dot tun- 
neling processes diagonal in the lead indices. The non- 
diagonal source-drain cotunneling may be treated by the 
canonical transformation method only perturbatively, 
but the adiabatic and non-adiabatic contributions into 
tunnel conductance may be sorted out also in this case. It 
will be shown that the time-dependent contribution into 
current-voltage characteristics of moving nano-object is 
significant near the boundaries of Coulomb diamonds. 
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II. MODEL 

We choose for the realization of ac-driven tunnel con- 
ductivity the simplest model of a nanoobject, which is 
widely used in the studies of single-electron tunneling 
(SET). A nanoobject is represented in this model by the 
quantum well with resonance level (see Fig. [U up- 
per panel). The SET regime arises due to the Coulomb 
blockade effect: addition energy for the second electron 
in a singly occupied dot is Sd + U 3> Tj , where Tj is the 
tunneling rate to the left (j = I) and right (j = r)lead, 
and U is the capacitive energy of the dot. To sepa- 
rate the time-dependent SET from the Kondo ZBA, we 
consider tunneling between spin-polarized leads, where 
spin-flip processes are inelastic, because the continuum 
of electron-hole pairs with the opposite spins in the leads 
responsible for the Kondo effect is gapped. Two types of 
spin polarized (magnetically ordered) metallic leads are 
presented schematically in Fig. [TJ In the middle panel 
the leads are formed by " half- metallic ferromagnets" 26 
with gapped spectrum of minority spin carries. In the 
lower panel characteristic for p-type degenerate dilute 
magnetic semiconductors^ the carriers are the minority 
spin holes. 

Virtual tunneling results in a shift of level positions 
in QD. This renormalization ("Friedel shift") is also spin 
dependent. As a result the spin polarization of QD ad- 
justs to that of the ferromagnetic lead (see calculations 
below). All Kondo processes are quenched in this regime. 

The Anderson Hamiltonian modeling SET has the 
form 

H = H A + H h + H tun , (2.1) 

where the terms 
ff d = 5> A |A)(A|, H h =J2J2 e jk^ jk , a jk, (2-2) 

A j=l*T ka 

describe the electron states in the isolated dot and two 
metallic (semiconductor) leads, respectively. We write 
Hd in terms of its eigenstates |A) (Hubbard representa- 
tion). This trick allows one to take all intradot interac- 
tions into account exactly even when the contact with 
the leads is switched on. The tunneling Hamiltonian 

H tna = J2(Vjk4a jk(T +K.c.). 

jka 

may be rewritten in the Hubbard representation by ex- 
panding the creation operator S a in terms of the config- 
uration change operators X AX = |A)(A| which connects 
the states in adjacent charge sectors with N and N — 1 
electrons in the dot. We confine ourselves with the sim- 
plest case, where only three charge sectors N = 0,1,2 
are involved in SET Hamiltonian. Then the Hamiltonian 




FIG. 1: Upper panel: model of quantum dot in time- 
dependent contact with leads. Middle and lower panels: den- 
sity of electron states in the leads, occupied and empty states 
in the dot for N = 1; mechanisms of spin-polarized electron 
cotunneling (middle panel) and hole cotunneling (lower panel) 
are indicated by arrows. Figures in parentheses point out the 
sequence of electron tunneling acts in cotunneling processes. 

(|2.1j) acquires the following form 

H = Sd J2x^ + E 2 X 22 + J2^a] k(T a jk!T 

<y jka 

+ £[^4^A 0CT + aA ff2 )+H.c.]. (2.3) 

jka 

Here the quantum numbers A, A = 0, er, 2 correspond to 
empty, singly and doubly occupied states of QD, respec- 
tively, E2 — 2sd + U is the energy of doubly occupied QD. 
The last term in this Hamiltonian is time-dependent. 



III. CANONICAL TRANSFORMATION OF 
ANDERSON HAMILTONIAN 

Our program is to exclude the tunneling term from the 
Hamiltonian (j2.3|) by means of the canonical transforma- 
tion 

H = e s He~ s , (3.1) 
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then derive the tunnel current operator in the new basis 
and calculate the tunnel conductance. It was shown in 
Ref. HU that this transformation may be performed ex- 
actly in the absence of Coulomb correlation, provided the 
energy level Ed falls into the energy gap and remains there 
after renormalization (Friedel shift). It will be shown be- 
low that the matrix S still may be found in the presence 
of Coulomb blockade under the same condition of dis- 
creteness of renormalized ci-level provided the spin-flip 
processes are quenched. As may be perceived from Fig. 
[TJ this condition is realized for the majority spin elec- 
trons in the case of half-metallic ferromagnet and for the 
minority spin holes in the case of p-type dilute magnetic 
semiconductor. 

As usual, the canonical transformation is made by 
means of the Bakcr-Hausdorff expansion 



'ffe- s = ]T-[S, [S,...[S,H]...]] 



(3.2) 



71 = 



In the non-interacting case, the second quantization op- 
erators in S and H possess Fermi-like commutation re- 
lations, the Hamiltonian H is a quadratic form and the 
tunnel operator conserves spin, so the series in the r.h.s. 
of (|3.2p may be summed exactly^ The Coulomb block- 
ade separates the Hilbert space for the dot electron oper- 
ators into charge sectors divided by the energy gaps. As 
a result these operators lose the simple Fermi statistics. 

We are interested in the strong Coulomb blockade 
case and start with the simplest case, where the ground 
state of the dot corresponds to N = 1 in the limit of 
U ^> \e,i — ep\, so that the doubly occupied states are 
completely suppressed. Then only the states A = er, 
are retained in the Hamiltonian (|2.3[) . In particular, the 
anticommutation relation for Hubbard operators mixing 
the adjacent sectors N = 0, 1 has the form 



A. Time-independent transformation 

In this section we generalize the canonical transforma- 
tion proposed for an Anderson model applied to tran- 
sition metal impurities in semiconductors ! 28 ' 29 In those 
calculations the intra-impurity Coulomb interaction was 
taken into account in Hartree approximation. Here we 
take the Coulomb blockade term exactly, when summing 
the series (|3.2I) . at least in the charge sectors N = 0, 1. 
The antihermitian operator S is looked for in the form 



jka 



u* jk X°°a 



jka 



If the spin flip processes are neglected, the canonical 
transformation is made for each spin projection sepa- 
rately. One may apply the transformation to the ma- 
jority spin in the case of electron tunneling and to the 
minority spin in the case of hole tunneling in two models 
introduced above. 

One easily derives the commutation relations 

jk 



[S a ,a] 



jka I 



~ U jka X 



rrO 



To shorten notations, we introduce the quantities C\ — 
u jkaa] ka and 7 2 = J2 u jkaU* kcr . Besides, we omit the 

jka jk 

lead index j and specify the band state by a single in- 
dex k characterizing both the lead and the electron wave 
number. Using these definitions and the above mentioned 
idempotency of operator K ai we obtain the following ex- 
pressions for the transformed operators 



X 



a0 



X a0 ,X 



Oct' 



X c 



x°%, 



(3.3) 



which follows from the obvious multiplication rule 
A AlAl A A2A2 = (5 Ai a 2 A AiA2 . Even disregarding the spin- 
flip processes each commutation operation in the ex- 
pansion (|3.2[) generates operators K a = X aa + X 00 . 
In spite of this, the Baker-Hausdorff series still can be 
summed exactly, because these operators are idempotent, 
K"K a — K 17 , like conventional fermion occupation op- 
erators d\d a . This summation will be performed in the 
following subsection. 

Time-dependent problem is more complicated because 
in that case the canonical transformation should be ap- 
plied to the operator 



H 



(3.4) 



In subsection B we will show that the canonical trans- 
formation is operational in this case as well, at least in 
some important special cases. 



oS-X^e' 8 ' = X a0 cos la + ClK a j- 1 sin la , (3.5) 



oL =e s "aLe s ° 



a\ a + u^CiK'j- 2 (cos 7 . - 1) - u^X* ^ 1 sin 7CT , 

(3.6) 

and 

e s Cle- s = Ct + C^A CT (cos7 a -l)-A CT0 7 CT sin 7ff , (3.7) 

The tunneling term in the transformed Hamiltonian is 
eliminated, provided 



Uko 



laV k 



with 



(ska - E d a) tan 7, 
dL a (e) 



= g ka V k (3.8) 



tan 7o- = — 



de 



(3.9) 



e—Eda 
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k £ £ka 

Then the transformed Hamiltonian takes the form 

H a = E da X°° + J] ~el k ,a\ a a kl(y (3.11) 

kk' 

where the renormalized level position is given by the 
equation 

P £ d + Z a tan 2 la - (T a + T*) tan 7g 

-C'dcr — — — 2 ' \ 6AZ ) 

1 + tan z 7 CT 

with Z ff = l~ 2 Y.k £ kU* krT u ka , and T CT = 7J 1 Efc(^fc"D- 
Substitution of Eqs. JO) - (j3~T0)l in (j3~T2)) transforms 
it into the conventional for m 28 i 29 

-Edtr = £d + L a (Eda)- (3.13) 

where the self energy L a (E) has only real part (Friedel 
shift), provided the level Ed„ remains within the gap, 
which is the case for a =] in a configuration of Fig. [TJ 
middle panel. 

A canonical transformation for the second spin com- 
ponent a =1 should be done more carefully, because the 
bare level £d falls into continuum of spin-down states. 
One should be accurate with turning to the thermody- 
namic limit, where the sum in the right-hand side of 



(13.10[) transforms into the integral and acquires the imag- 
inary part thus making the Hamiltonian non-Hermitian. 
The recipe is to keep the spectrum of electrons in the 
leads discreet when doing the canonical transformation. 
Then equation (|3 . 13[) has JV + 1 solutions e$, where Af is 
the number of state in the valence band e k ^. Using Eqs. 
(|3.10p and (|3.9|) . the corresponding coefficients -fi may be 
found. In accordance with pil A[) . the factor cos7i deter- 
mines the weight of the (i-component of hybridized wave 
function of dot electron in the state i. One may identify 
the state i = m producing the maximum value of cos7i 
with the center of future 'Friedel resonance', which arises 
in the thermodynamic limit J\f — > oo. In this sense the 
state Edi = £mi is formally defined from the equation 
E m | = Ed + Li(e m i). This level is shown by the dashed 
line in the middle panel of Fig. [TJ It corresponds to 
the bunch of excited states of QD coupled to the leads 
with N = 1 and spin oriented antiparallel to that of the 
magnetized leads. 

The spectrum of continuous part of the Hamiltonian 
(|3.11[) is determined by the expression 

el k , = e ka S kk , + W£ k , (3.14) 

containing the "scattering" matrix element which even- 
tually predetermines the tunnel current. 
The general form of this matrix element is 



W£ k , = 



u kaU k , a 



'jL K a 



(2e d - Ska - - COS7 CT ) - 



T a + T* 2 cos 3 7g - 3 cos 2 7<T + 1" 
2 cos 7cr sin 7^ 



K°{V k ui, a + u ka V k *,) S -^. 

7<r 



(3.15) 



Using Eqs. JSHJ and (l3~T2"j) we get 

1 



T a = t; 



tan7 CT 



L a {Eda) 



and after some algebraic manipulations the scattering 
matrix element is eventually transformed into a quite 
compact expression 

wl k ,= 



V k V k *,K° 
where 



A fc „ A fc 



1 "N J3l \ ^^(Edc) 2/ s 
R M + A A R (7<t) 



A ka A k , 



1/1 + tan 2 7(T - 1 



tan 2 7c 



(3.16) 
(3.17) 



Afc ff = Eda — £fc<r- Equation pil A| holds in the static 
case or, as we will see below, for adiabatically slow time 
variations of the tunneling amplitudes. If we want to 
study nonadiabatic corrections a more general formula 
Eq. pil A|) should be used. 



Equation (jlll A[) generalizes the familiar 2-nd order 
expression for the single electron tunneling amplitude 
through the QD, which takes into account both renormal- 
ization of the energy level of dot electron (|3.13p and re- 
construction of the band continuum (jlll Ap . Far from the 
resonance tunneling regime (in the center of the Coulomb 
diamond diagram, see e.g. the middle panel of Fig. [JJ 
one may neglect the second term in the square bracket of 
Eq. pil Ap . and the tunneling matrix element acquires 
the simple form W° k rk , = Ji k ,rk' K a i?( 7(T ), where the 
first factor is the off-diagonal matrix element of indirect 
exchange between the leads and the dot due to electron 
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cotunneling, the second and the third factors regulate 
the occupation of the dot level and the normalization of 
electron wave functions, respectively. The transparency 
of QD is, of course, exponentially weak 

Now we turn to the resonance regime illustrated by Fig. 
[21 where the level e d — ed + U driven by the gate voltage 
v g approaches the level ep from above. The level e\ may 




t e 



N=l N=2 




drain 



FIG. 2: Cotunneling mechanism in the resonance regime for 
second electron in QD on the boundary of Coulomb window. 

be occupied only by a spin-down electron. Since it falls 
into the gap of spin-down density of states in the leads, 
the canonical transformation introduced above may be 
performed in a similar way, provided the addition energy 
for the first electron sa e d falls deep enough below €f 
and the corresponding processes are suppressed. In this 
regime A, A = a, 2, and one should use the commutation 
relations 



[X a2 ,X 2a '} + = x aa ' 



(3.18) 



instead of (|3.3p . Correspondingly, one should insert 

K a _ X aa + X 22 m equat j ons 



X 2S = aX 2s cos-/ a 



7<x Smjcr, 



(3.19) 



substituting for Sd and the correlation function K a taken 
from (|3.18|) . It should be stressed that the tunneling 
through the QD is impossible at zero source-drain bias 
because of the spin blockade in spin-polarized electrodes. 
Only spin-up carriers exist around Fermi level, and these 
electrons may be injected into QD only when accompa- 
nied by the spin-flip excitations given by the operators 
X** in the intermediate state with TV — 1 of cotunneling 
process. These processes are inelastic and exponentially 
weak (~ V 2 in transparency and ~ V A in conductance). 
More detailed discussion of spin-dependent tunneling is 
postponed till Section [TV Al 



B. Time-Dependent Transformation 

As was discussed above an experimental conditions 
may be created when the tunneling matrix element Vk (t) 
in the Hamiltonian (|2.3p becomes time dependent. The 
canonical transformation (|3.ip . as described in the pre- 
vious section, cannot be straightforwardly applied. Its 
generalization is in order. 

We start with the temporal Schroedinger equation 



Ctp = 



(3.22) 



and look for the time dependent transformation matrix 
S(t), which transforms (|3.22p into 



Ctp = 



(3.23) 



with transformed operator 

C{t) = e s «tfe- s « - ihe s ^^e- s M 

for the new wave function ip = e s ^^p. 

The Hamiltonian H satisfying the equation (|3 . 23[) can 
be written as 



H = e s He- s + ih d\e xs Se- xs . 



(3.24) 



«L + u* ka CtK^- 2 (cos la - 1) - ul.aX 2 ^- 1 sin 7 , 

(3.20) 

for transformed creation operators. 

Then the transformed Hamiltonian for spin-down elec- 
trons has the form 



H i =E dl X^+J2^kk'4l a k'l 

kk' 



(3.21) 



where E d i and e1 k , are given by the same Eqs. (|3.13p and 
p,14p as in the previous case, but with the energy 



(see Appendix). 

The Hamiltonian (|3.24p contains now two terms of 
which the first one is just a modification of the Hamilto- 
nian (|3.ip of the time independent case. It means that 
all the equations in Subsection IIII Al hold except for Eq. 
(|3.8p . which defines the coefficients Uka of the canonical 
transformation. These coefficients must be now found 
anew. There is also the second term in the right-hand 
side of Eq. (|3.24p . which is responsible for non-adiabatic 
effects. 

In order to find the canonical transformation parame- 
ters Uku we write explicitly the condition 
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Ukala 1 l( £ d - Z a ) sin T CT cos 7 CT - (e k<7 - Z„) sin j a + T a cos 2 7 CT - T* sin 2 7 - T a cos 7 CT ] + V fc cos 7 CT = 



fe<7 



sinTo- 

7<T 



E 

fc' 



u k > a u k , 



■(sin7 CT cos7, T +7,, 



2sin7 CT ) + u* k , a u kla — (1 
27 2 



Sill 7cr COS 7cr . 

7<r 



(3.25) 



of elimination of the QD - lead tunneling in the trans- 
formed Hamiltonian. Here both the tunneling amplitude 
Vfecr and transformation parameter u ka are functions of 
time. The condition pil Bp contains a number of terms 
with the time derivatives ii k(T . Neglecting these time 
derivatives would correspond to the adiabatic approxi- 
mation where the variation of the tunneling amplitude is 
slow enough and the whole electron system always have 
enough time to readjust to the varying tunneling ampli- 
tude without additional level mixing. Then one can check 
straightforwardly that Eq. (|3.8|) with the time dependent 
Vfe(t) solves Eq. pTlBl) . 

Now we carry out a more general analysis going beyond 
the adiabatic approximation. For this sake we multiply 



Eq. (|fflB|) by u 

equation 



ka 



and sum over k, which leads to the 



(e<j - Z a ) tan -y a + T a - T* tan 2 7^ = 



-ih^tan^-ih ^^ ^(u k , cr u%, a -u* k , tT u k ' t r)- (3.26) 



Substituting Eq. B|) into Eq. pil Bp we get the equa- 
tion 



u ka [{e k a - Z a ) tan7 CT + T a ] - 7V* 



ka 



tan 7 CT 



-u ka 



k' 



l k'a'a u k', 



+ ihu kcr tan7 CT . (3.27) 



Neglecting the time derivatives of transformation param- 
eter u k (t) in the r.h.s of Eq. (jlll Bp . which corresponds 
to the adiabatic approximation, yields Eq. Ij3.8jl with 
the time dependent tunneling amplitude Vk(t). All the 
other equations obtained in Subsection IIII Al also hold. 
Accounting for the r.h.s. of Eq. pil Bp allows one to 
obtain nonadiabatic corrections. 

Having in mind that the quantities e d , 7<t and Z a are 
explicitly real, we separate real and imaginary parts in 
Eq. piIBp and thus get two equations 



and 



T a - t; = -2ihi« 



rp I rp* 

(s d - Z a ) tan 7, + " a (1 - tan 2 la ) 



(3.28) 



tan7 CT ^ „ 

~ o > [Uk'aUk'a-U klcr Uk'a)- 



27 2 



k' 



(3.29) 



These equations may be instrumental in looking for solu- 
tions u krT (t) for specific problems. It follows from (|3.28[) 
that the quantity T which was real for the time indepen- 
dent case, remains real also in the adiabatic approxima- 
tion for the time dependent case. 

Returning back to the transformed Hamiltonian (|3.24D , 
we note that the first term e s He~ s is now time- 
dependent due to the time dependence of S. Carrying 
out the transformation in the same fashion as in the pre- 
vious section we get the time dependent energy level 



E. 



(a) 
da 



E. 



(b) 
da 



where 



p (o) _ £ d + Z a tan 7 CT 

&da " 



(Ta +T*)tan7 CT 



(3.30) 



(3.31) 



1 + tan^ 7(t 

is the same energy level (|3 . 1 2|) as before but with the coef- 
ficients Z„ , T a , T* , 7 CT depending parametrically on time 
t. Thus, in adiabatic approximation the time dependence 
of the resonance level position E d a J (t) is determined by 
Eq. I3.13P with the time-dependent self energy part 

\v k (t)\ 2 



L a (e,t) 



k 



Ska 



(3.32) 



There is also the nonadiabatic correction 

2 



E 



(<) = --lh — 2_JU ka aU k 

~ ka 



Mm - 

da 



u kaa u kaa) 



(3.33) 

To calculate the time dependent coefficients, one has to 
specify the form of tunneling amplitudes V k (t). An ex- 
ample of time-dependent tunneling will be considered in 
the next section. 



IV. TUNNEL CONDUCTANCE 

In this section we study the tunnel conductance basing 
on the transformed Hamiltonian H. The tunnel current 
may be calculated, e.g., by means of the Keldysh tech- 
nique, where the bias eV is included in the zero order 
Hamiltonian and the scattering ~ Wik, r k' is considered 
as a perturbation. First, we calculate the spin-polarized 
current through an immovable quantum dot and then 
discuss the modulation of this current due to oscillatory 
motion of the dot. 
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A. Tunneling through static dot 

Far from the boundary between two adjacent charge 
sectors with N = 1 and N = 2, where both the levels 



and e d 2 ^ are far from the chemical potential fi, the 
Keldysh method applied to the Hamiltonian (|3.21| in a 
single loop approximation gives the conventional golden 
rule equation 



J 



I = e Y^2^2 \ W ilrk> \ 2 S(eika +ev- e rkr<T ){f(e lka )[l - /(e rfc , CT )] 

kk' a 



f(erk>*)[l-f(eik*)}} 



(4.1) 



r 



where ev is the source-drain bias, f(tko) is the Fermi 
distribution function for spin-polarized electrons, and the 
scattering amplitude Wik, r k' is defined in Eq. pil A[) . 

The standard Coulomb diamond diagram for tunneling 
conductance G(v g ,ev) is distorted in the region of gate 
voltages v g corresponding to the change of QD occupa- 
tion (N = 1) — > (N = 2). In case of completely spin 
polarized dots the Coulomb step corresponding to the 
resonance E(N = 1, v g ) = E(N — 2, v g ) — fx in the cur- 
rent voltage characteristics is absent at zero bias and zero 
temperature due to the spin blockade. The chemical po- 
tential fi is pinned to the Fermi level of spin up electrons 
eft, but the resonance level E di = E(N = 2)-E(N = 1) 
belongs to the down spin electron (see Fig. Thus the 
tunneling at zero bias ev = is suppressed by the spin 
blockade. This blockade may be surmounted by means 
of finite source-drain bias compensating the energy gap. 
However, the conditions for the onset of spin up and spin 
down tunnel current are different. 



N=2 




N=l 




ev 



FIG. 3: Coulomb diamonds diagram for tunneling conduc- 
tance G(v g ,ev) near the border (N — 1)/(N — 2) for down- 
spin and up-spin electrons (solid and dashed lines, respec- 
tively) . Three bars on the v g axis mark the values of the gate 
voltage corresponding to the level Edi in resonance with e c , 
ef, £v (bottom- up) . 

The boundary of the Coulomb diamond with N = 1 
for spin down electrons follows the evolution of the res- 
onance level Edi un til this level approaches from above 
the bottom of conduction band e c i . When it crosses the 
band edge, the resonance tunneling is no more possible at 
eV < A CVl so the Coulomb resonance line for G(v g ,ev) 



deviates from the linear behavior, as it is sketched in 
Fig. [3] (solid line). Negative-bias part of the Coulomb 
diamond corresponds to the hole tunneling through the 
occupied resonance level in the sector N = 2. It is dis- 
torted in the same way in the region of v g where this level 
matches the top of the valence band (second branch of 
the solid line). When the down-spin level is deep enough 
in the conduction or valence band, the linear behavior is 
restored again. 

The blockade for spin up electrons near the boundary 
(N = 1)/(A = 2) is lifted at ev compensating the spin 
flip excitation in the dot. As was mentioned above, the 
spin up electron tunneling is allowed provided the dot 
is excited to the state Edi gi yen by the solution of Eq. 
(|3.13p for the doubly occupied dot, and the splitting en- 
ergy may be estimated as Afj, « \L{Ed]) — L{Ed{)\ (see 
Ref. [25| for experimental determination of such splitting 
in the Kondo tunneling regime). The line of resonance 
tunneling for spin-up electrons is drawn in Fig. [3] by the 
dashed curve. 

Besides, the tunneling transparency is especially sensi- 
tive to the position of the dot level in the near vicinity of 
band edges. To investigate this dependence, let us find 
the explicit equation for the tunnel conductance from 
Eqs. (|4.1|) and (|III A[) . Changing summation over kk' 
for integration over ee' in a usual way and performing 
standard calculations, one gets the equation 



G l (v g ,A cv 



T,T r 



h (2tt 2 ) \Edi - e c 



R(jl)L(Edi) 



(4.2) 



for the threshold value of ev — > A cv +0 in case of 2D elec- 
tron gas in the planar leads. Here Ti = is the 
tunneling rate for the lead i obtained in the approxima- 
tion of Vik = Vi and constant electron density of states 
Si(e) = Si, which is valid for 2D electrons. It should 
be taken into account that the definition of the posi- 
tion of the resonance level Edi falling into the conduction 
band continuum implies the procedure described below 
Eq. (|3.13p The function L(e) has a logarithmic singular- 
ity at the band edge, ReL(e — > e c ) ~ — ln(|e — e c |/A c „). 
As a result, Eq. Q3.13P has either one or two solutions 
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-(2) 



v g relative to 



depending on a position of the level e d < 
the band edge£L 

The resonance factor \Edi — £ c \ 2 in the denomina- 
tor and singular factor L(Edi) in the numerator of Eq. 
(|IV A[) result in noticeable enhancement of G at the 
boundary of Coulomb diamond. 

This enhancement is characterized by evolution of the 
ratio 



p{v g 



\E 



(4.3) 



from its value in the middle of Coulomb diamond to that 
in the vicinity the point Ed{v g ) — e c . Here and below we 
omit the superscript (2) in the notation of Ed- Far from 
the Coulomb resonance the difference between \Ed\\ and 
Ed is small and tan 2 7j ~ V 2 /(£d — £ c ) <C 1. The function 
R(jl) tends to 1/2 in this limit, so that the factor p may 
be estimated as p « 1/4(6^ — e c ). Near the band edge 
e c the factor i?(7j.) « cot7| = \/\E d i — E c \/T r [see Eqs. 
(i331) . (j3~T71) ]. so that p ~ [\E d[ ~ e c \T r }- 1 . Numerical 
estimates of this enhancement are presented at Fig. 2) 
In these calculations the density of states was assumed 
to be constant in the lower part of 2D conduction band, 
so that the self energy in the r.h.s. of Eq. (|3. 1 3[) may 
be approximated as L(E) ss T r Ln(E/D), where D is the 
effective width of conduction band, and the argument E 
is complex. Then the solutions E — Ed[ of Eq. (|3.13|) 
are expressed via the Lambert W- function W(n, x) 



E, 



di 



-T r W(n,-r- 1 e- ed/rr ) 



(4.4) 



(here the reference point is e c — 0, all energy parameters 
are measured in units D, index n enumerates branches 
of W- function) . There are two solutions for Edi near the 
band edge^i The principal branch n = gives the dis- 
crete level in the gap, and the branch n = 1 corresponds 
to the resonance in the band. Here we are interested in 
the latter state. 

The amplification reaches its maximum when the dif- 
ference \Edi —e c \ comes up with T r , therefore the smaller 
r r , the bigger is the enhancement factor p. The factor 
RL I A cv in the square brackets in Eq. (jIV A[) behaves as 
e 1/,2 lne in the vicinity e = \Ed{ — £ c \ of the band edge 
and thus gives additional contribution to this enhance- 
ment. Both analytical and numerical estimates confirm 
this statement. Experimentally, this effect should be ob- 
served as an increase of tunneling conductance on the 
boundary of Coulomb diamond in the vicinity of the 
threshold value of eu t h = A cv . A similar effect should 
arise on the hole side of the Coulomb diamond diagram 
(et>th = —Acv), where the occupied level Ed\_ crosses the 
top of the down spin valence band (solid lines in Fig. [3]). 



B. Tunneling through moving dot 

As is shown above, the canonical transformation al- 
lows one to distinguish between the adiabatic and non- 
adiabatic contributions to the tunneling amplitude. Let 




0.0 0.1 0.2 

Prime dot level 



FIG. 4: Amplification factor p(e^ ) for two values of tunnel- 
ing rate T/D = 0.01 (lower curve) and F/D = 0.02 (upper 
curve) . 



us first discuss the adiabatic corrections to the inelastic 
current given by Eq. (|4.1|) and illustrated by Fig. 
In the vicinity in the point (N — l)/(N — 2) of the 
Coulomb diamond diagram the adiabatic position of the 
down spin level E d a ) given by the solution of equation 



E d f(t) = 



E d\ (*) _ e lk 



(4.5) 



rocks around e cr or e v i (see Fig. [2]). This solution de- 
pends parametrically on time via the oscillating tun- 
nel coupling |Vjfc(£)| 2 . The above analysis of the static 
case prompts that time-dependent corrections become 
significant at finite bias close to the threshold value 

(a) 



±A C „. The adiabatic evolution of the level E' d ^ 



e^th 

in time near e c is given by the same Eq. (]4.4|) . where 
the tunneling rate T r parametrically depends on t. If the 
time-dependent perturbation is weak in comparison with 
the static value of tunneling rate, the temporal compo- 
nent of E d ? may be found perturbatively. Representing 
the tunneling rate as T r = T r Q + ST(t) and expanding 
Eq. (|4.4p around the time-independent value marked by 
index '0', we get 



E d a ?[t) = E i 



(o) 



Wo 



1 + Wn 



Ed 



+ W 5T(t) (4.6) 



When deriving Eq. (14. 6|) . the equality W{x) — x[l + 
W{x)]dW/dx is used. 

This time dependence turns into the corresponding 
adiabatic time dependence of tunnel conductance, mainly 
via the enhancement factor p(Edi) (|4.3[) . One may ex- 
pect that the slow adiabatic variations of G(t) will be 
especially distinct at v g corresponding to steep slopes 
of p(t) (Fig. 0|) at ev ~ ±A CU near the boundary 
(N = 1)/(N = 2) of the Coulomb diamond diagram 
(Fig.©. 
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C. Weak time dependent perturbation 

Calculation accounting for nonadiabatic corrections to 
tunnel conductance is generally an extremely compli- 
cated nonlinear problem. To make it tractable, we as- 
sume that the time dependent part of VjUcr is only a 
small periodic perturbation with respect to the time inde- 
pendent part Vj^l and consider only the linear response 
given by the first harmonics. The nonlinear effects will 
be discussed separately This approach allows one to pick 
out first nonadiabatic corrections to tunneling amplitude 
W kk > i which turns out to be small everywhere except in 
the vicinity of band edges. 

Let us assume that the tunneling integral has the form 



Vk(t) = V h 



- T/(0) 



(1) 



COS Lot 



(4.7) 



(here and below indices ja are omitted for the sake of 



brevity). Here V fc (0) > V^' . The solution of Eq. (ImB) 



r(l) 



is looked for in the form 



Uk 



,(o) 



(!) + 
u k COS Lot 



sin ujt 



(4.8) 



and collect separately all the terms containing cos Lot and 
sin cot respectively. Substituting then u k in A|) . we 
reduce the scattering amplitude W kk i to the following 
form 



W kk , = W 



kk' 



W^cosLot + hioW^, sin cot (4.9) 



where the coefficients of the cosine and sine terms can be 
explicitly calculated. The cosine coefficient 



SW kk r (1) 

sv q q 



sw kk 



sv* q 



(4.10) 



is obtained by varying equation A[) over the hy- 
bridization potential V q , whereas the sine coefficient 



5W kk 



sw kk 



-W, 



(2) ,nonad 
kk' 



(4.11) 



where the time dependent corrections to are also 
small. Then we vary all the coefficients in Eq. (|3.28|) 
with respect to = g k V^ 0) [see Eq. ([3T5]) ] and 



is obtained by varying equation (|III A[) over the transfor- 
mation parameter u q . After the variation u k in the form 
Eq. p.8p can be substituted. Then 



|^(2) ,nonad 



v^v k {0) 



v fc (1 V fc ( , 0) * 



A ka A k , 



;/ (0) T/ (0)* (1) (0)* T/ (o)* T/ (o) 

V k V k> n2(~, \ \ " V k" V k" ~~ V k" V k" 



r 



with 



Rn 



Rn 



i/T+ tan 2 -f a 



The most divergent terms in (|4.10[) and (jIV C[) appear 
when we vary only explicitly written in Eq. ((Ill A[) 
and uy in Eq. (jlll A|) . As a result keeping the leading 
terms we have 



The corresponding corrections to the tunneling trans- 
parency near the conductance band edge behave as 



SGf\v g ,A c 



e 2 r<r« 



R 2 (ll) 



and 



h (2^ 2 ) \E dL - e c 



SGf'(v g> A cv ) 



cos cot (4.14) 



1 1 



A ka A k ' fj 



(4.12) 



and 



(0) T/ (1)* t/(0)* t/ (1)> 



(hoj) 



e 2 r,if } i? 2 ( 7i ) 

h (2^ 2 ) - £c p/2 2Vr, 



: sincji 



where 



(4.15) 



is the adiabatic correction to the tunneling rate for m = 1 
and the correction due to weak non-adiabatic effect for 
m = 2. Equation (|IVC|) uses the fact that 



RiTfa 



A L A 2 V ; tan 7(T 



(4.13) 



71 



and tan 7^ 



\Edi - £c 
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for Edi — > e c when \E dl - 

r{m) ^ y(m) 



: c \ ?s r r . We use here the 
similarly to the one used in 



approximation V fe 
Eq. (ITVAl . 

The term with the sine in equation (|4.9[) causes a phase 
shift between the oscillations of the hybridization param- 
eter and the resulting current through the dot. If we 
neglect the k dependence of the coefficients (14. 10)) and 
(ITVCj) this phase shift can be readily found 



■SGf\v g ,A cv ) 



tp = — arctan 



It is expected to be generally rather small due to the fac- 
tor gkfiuj, which is usually very small unless the dot level 
approaches the band edge e c . However close to the edge 
this factor diverges and results in an increasing phase 
shift, 



(p w — arctan ■ 



\/^r\Edi — £c 



— arctan 



As result the phase shift ip may become essential at 
Edi — * e c - 

The non-adiabatic corrections to tunneling conduc- 
tance acquire the simple form of phase shift in oscillating 
cosine function only until the parameter fk — ghh^ ~ 
Hcu/A is small and the perturbative approach is valid. 
With increasing / one may expect appearance of higher 
harmonics ntv in oscillating conductance. With further 
increase of the amplitude V 1 - 1 ^ the language of quasi- 
energy levels 3 -^ is more appropriate. We plan to discuss 
it elsewhere. 



V. CONCLUDING REMARKS 

We have discussed in this paper a new approach to the 
Anderson model for a half-metallic electron liquid in a 
tunnel contact with a moving nano-shuttle under strong 
Coulomb blockade. It is shown that in the situation 
where the spin-flip cotunncling processes are suppressed 
at low energies, the exact canonical transformation elim- 
inating the tunneling term in the Anderson Hamiltonian 
exists even in the presence of strong Hubbard repulsion 
in the shuttle and time-dependent lead-shuttle tunnel- 
ing. This canonical transformation in principle allows 
one to sort out the slow adiabatic renormalization of the 
energy levels and the tunnel transparency and to con- 
sider non-adiabatic corrections at least perturbatively. 
One may also include the inelastic spin-flip processes in 
the canonical transformation in the 4-th order of pertur- 
bation theory in Vi, but these weak corrections do not 
change the above qualitative picture. We also have cal- 
culated the weak non-adiabatic corrections to tunneling 
transparency in a specific model of periodic sinusoidal 
motion of the shuttle. Undoubtedly, strong non-adiabatic 
effects should be taken into account in a more refined 



scheme: in the case of a periodic time-dependent pertur- 
bation one should appeal to the Floquet theorem in the 
time domain and use the quasi energy language. 32 ' 33 

Because of the energy gap for spin-flip processes in 
the half-metallic leads the Kondo cotunncling processes 
are suppressed, and the zero-bias anomaly in the tunnel 
conductance is absent. Physical manifestations of the 
shuttling mechanism under discussion arise in the form 
of time-dependent enhancement of conductance at finite 
bias near the boundaries of Coulomb diamonds on the 
phase diagram G(v g , ev). The boundaries themselves are 
distorted due to the complete spin polarization of carriers 
(see Fig. [3]), and even the Coulomb blockade step corre- 
sponding to occupation change from odd to even number 
of electrons in the dot is absent at zero bias. One may 
expect appearance of quasi energy satellites in this part 
of the phase diagram G(v g , ev) when the shuttle motion 
is essentially non-adiabatic. This regime is a subject for 
future studies. 



APPENDIX A: TIME DEPENDENT 
CANONICAL TRANSFORMATION 



To diagonalize the Schrodinger equation (|3 . 23[) for the 
wave function ijj = e s ^ip : we conjecture the form 



H = e s He~ s + S 1 



for the canonically transformed Hamiltonian. Then sub- 
stituting H and i\> into Eq. (|3.23|) . we obtain 



d 



dt 



(e b He~ b +S!)e b ip . (Al) 



The time derivative of exponent is found by means of the 
operator equation 



-re" 
dt 



e ^s e -Hs-i) dX = e s 



e - xs Se xs dX 



(A2) 



Then using (fM| and (IX2^|) in (fM|) yields 

e s V + e S Hip = e s Hijj + S x e s ip. 



ih 



dXe xs Se~ xs 



U=o 



from where we straightforwardly get Si and Eq. (|3.24D . 



APPENDIX B: TUNNELING AMPLITUDE FOR 
WEAK PERIODIC POTENTIAL 



To derive the first non-vanishing time-dependent cor- 
rections to tunneling procedure, we vary the terms in Eq. 
(|3.28p . This variation procedure results in the two linear 
equations 
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(i)V(o) (i)t/(o)* fc ft A°)* (!) (i)*t/(o)"\ fiA°)*iA 1 ) 

I 



(Bl) 



where is explicitly real as defined in Eq. (|37 

The first of these equations contains only explicitly 
imaginary terms whereas the second one contains only 
real terms. Therefore we have only two equation for four 
real variables (two complex variables). It is readily veri- 
fied that 



4 1} 



9k 



1 - (hug k y 
1 - {hug k ) 2 



(i) 



V, 



(1) 



(B2) 



is a particular solution of the set of equations (|B1[) . One 
can also see that in the limit g k hio — > this solution re- 
duces to the adiabatic approximation where we can take 
the results obtained in Subsection IIII Al for the time in- 
dependent problem and substitute there the tunneling 
amplitude ()4.7j) slowly varying in time. 



The general solution of the corresponding homogenous 
set of equations reads 



fuvg k 



hug k 



(Cife -I- iC 2k ) + [ hu>g k - 



^9k) 2V, 



Vr 



(0) 



,(0) 



-(Cifc — iCik)-, 



(B3) 



„(i) 



C lk + iG 



2 k 



It should be added to the particular solution (|B2[1 . 
The parameters C\ k and C 2k are meanwhile arbitrary. 
These parameters determined from the original cance- 
lation equations ((III Bp - piIB[) turn out to me small, 
C\, C 2 ~ (g k hu>) 3 and may result at most in nonadia- 
batic corrections ~ (g k huj) 2 . 

In order to calculate the tunnel current we need the 
scattering matrix element 

W kk . = W$ + W^ ad (B4) 

where the first term W kk , is obtained by substituting 
the solution P~5j) . (lB2|) into Eq. pil Ap and keeping the 
terms linear in g k hiu. This contribution remains finite in 



the adiabatic limit g k htu — > 0. However it contains also 
nonadiabatic corrections due to the third term in (|4.8p . 
These corrections are determined by the equation 

W^P ad = ih{u%,u k - u k ut,) 1 ~ C 2 ° S7 (B5) 

T 

, (COS7 - l) 2 v^, „ , 

+ lflU k U k ^_Ju k nU k „ -U k ,,U k n). 

~ fc" 

which contains derivatives and, hence, is purely nona- 
diabatic, i.e. disappears in the limit g k huj — > 0. Keeping 
only terms linear in g k huj in the scattering matrix ele- 
ment p34p . we come to Eq. (|4.9[) . 
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